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Abstract. We study properties of a single- component (spin polarized) degenerate 
dipolar Fermi gas of ^^^Dy atoms using a hydrodynamic description. Under 
axially-symmetric trapping wc suggest reduced one- (ID) and two-dimensional (2D) 
description of the same for cigar and disk shapes, respectively. In addition to a 
complete numerical solution of the hydrodynamic model we also consider a variational 
approximation of the same. For a trapped system under appropriate conditions, the 
variational approximation as well as the reduced ID and 2D models are found to yield 
results for shape, size and chemical potential of the system in agreement with the 
full numerical solution of the three-dimensional (3D) model. For the uniform system 
we consider anisotropic sound propagation in 3D. An analytical result for anisotropic 
sound propagation in uniform dipolar degenerate Fermi gas is found to be in agreement 
with results of numerical simulation in 3D. 



PACS numbers: 67.85.Lm,03.75.Ss,05.30.Fk,67.10.Db 



Study of a degenerate dipolar Fermi gas of ^^^Dy atoms 



2 



1. Introduction 

After the experimental realization of Bose-Einstein condensate (BEC) of ^^Cr [U |2], 
^^^Dy ||3j, and ^^^Er |1] atoms with large magnetic dipolar interaction, there have been 
renewed interest in the study of cold atoms, both theoretically and experimentally. The 
anisotropic long-range dipolar interaction acting in all partial waves in these atoms is 
basically different from isotropic short-range S-wave interaction acting in nondipolar 
atoms. Polar bosonic molecules with much larger permanent electric dipole moment 
are also candidates for possible BEC experiments |5]. Dipolar BECs have novel distinct 
properties. Due to the anisotropic dipolar interaction, the stability of a dipolar BEC 
depends on the scattering length as well as the trap geometry [2 |6l [7j. The shock 
and sound waves also propagate anisotropically in a dipolar BEC [HI E]- Anisotropic 
collapse has been observed and studied in a dipolar BEC of ^^Cr atoms ^0\. Anisotropic 
rotons [11] and anisotropic critical superfluid velocity [12] have been suggested and 
studied in dipolar BECs. Distinct stable checkerboard, stripe, and star configurations 
in dipolar BECs have been identified in a two-dimensional (2D) optical lattice as stable 
Mott insulator [13] as well as superfluid soliton [H] states. Anisotropic solitons in 2D 
have also been suggested in dipolar BECs [15]. A new possibility of studying universal 
properties of dipolar BECs at unitarity has been suggested \16\ . 

After the realization of BEC of alkali-metal atoms, degenerate nondipolar gas of 
fermionic ^Li [17], [18], and ^^Sr [19] atoms were observed. Later, superfluid states 
of paired ^Li [20] and [2T] Fermi atoms have been investigated. More recently a 
degenerate dipolar gas of fermionic ^^^Dy atoms with large magnetic dipole moment 
has been created and studied [22]- Realization of quantum degeneracy in ^^^Dy atoms 
should be considered as a doorway for the study of anisotropic superfluidity in dipolar 
fermions. Fermionic polar molecules, such as ^°K®^Rb, of large permanent electric dipole 
moment are also considered for this purpose [23]. The ^°K^''Rb molecule in the singlet 
rovibrational ground state has an electric dipole moment of 0.6 Debye, thus leading 
to a dipolar interaction larger than in the case of ^^^Dy atoms by more than an order 
of magnitude [1]. One advantage of studying the effect of dipolar interaction in a 
degenerate dipolar Fermi gas over that in a dipolar BEC is the remarkable stability of the 
degenerate Fermi gas. A BEC is usually fragile or short-lived for many experiments due 
to three-body loss through molecule formation. The three-body loss is highly suppressed 
in a degenerate Fermi gas due to Pauli repulsion among identical fermions. On the other 
hand a theoretical study of a BEC is much simpler than that of a degenerate Fermi gas 

due to the existence of a simple order parameter and a simple mean-field Gross- 
Pitaevskii equation in the former case. 

A microscopic description of a degenerate Fermi gas is complicated due to the 
difficulties with the antisymmetrization of a many-fermion system and finding an 
appropriate simple order parameter. Drastic approximations are often necessary to 
achieve this goal for fermions. There have been several different theoretical descriptions 
for a degenerate nondipolar Fermi gas [25]. There have also been a few studies of 
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the degenerate dipolar Fermi gas employing different types of approximations [26]. 
However, for a description of some macroscopic observables of a degenerate Fermi 
gas, a simple hydrodynamic description [27j often seems appropriate which does not 
require a precise antisymmetrization of the dynamics [211 I2H1 [29]. The effect of 
antisymmetrization is included approximately via an appropriate energy or Lagrangian 
density. A minimization of the energy leads to a well-founded variational approximation 
|3Uj . An Euler-Lagrange equation derived from such a Lagrangian provides further 
improvement over the variational results. Here we present a description of the degenerate 
spin polarized Fermi gas based on such a hydrodynamic model [271 128] . 

The present theoretical formulation for the degenerate dipolar Fermi gas starts 
in section 2.1 with the standard equations of classical hydrodynamics [27] with the 
appropriate equation of state including the kinetic energy of fermions filling the Fermi 
sea and the dipolar interaction among them. The equivalent of the Thomas- Fermi (TF) 
approximation for the dipolar gas is obtained by setting the velocity field equal to zero in 
the hydrodynamic equations. An equivalent classical energy density is written using the 
local-density approximation (LDA) [23]. A quantum pressure term is then introduced in 
the LDA energy density. With such a quantum pressure term, a nonlinear Schrodinger- 
type equation is derived for the density of fermions. For a moderate number of fermions 
(greater than 100 or so), the quantum pressure term is negligibly small and hence has 
insignificant effect on the result. However, the inclusion of the quantum pressure 
term allows one to write a dynamical equation for fermions in the form of a nonlinear 
Schrodinger equation to study the dynamics. The LDA or the TF approximation, 
on the other hand, allows only the investigation of static properties of the system. 
For example, these approximations cannot be used to study the sound propagation 
dynamics in fermions as in section 3.2 In section 2.2, we present a Gaussian variational 



approximation for the problem described by the LDA energy density. In section 2.3 



simplified models are derived in reduced dimensions, appropriate for cigar and disk 
shapes of the degenerate dipolar gas when there is a strong trap in radial or in axial 
directions, respectively. In section 2.4, using the hydrodynamic model, we obtain an 
analytic expression for the anisotropic sound velocity in the degenerate dipolar gas. 
In section 3.1, we present numerical results for stationary properties — shape, size, 
and chemical potential — of a trapped degenerate Fermi gas of ^^^Dy atoms, and 
compare with results from appropriate models in reduced dimensions and variational 
approximation. In section 3^ the anisotropic sound propagation in an infinite dipolar 
degenerate gas is demonstrated numerically and the velocities so obtained are compared 
with the analytical results of section [2] D. Finally, in section IV we present a summary 
of our study. 
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2. Analytical Consideration 
2.1. Hydrodynamic Model 

The normal one-component Fermi gas can be in the coUisionless regime where collision 
is rare or in the coUisional hydrodynamic regime where frequent collision due to the 
dipolar interaction allow the system to settle to local equilibrium, where the system can 
be described by simple hydrodynamic equation rather than a detailed multi-particle 
description. Here we consider the system in such a configuration. A semi-quantitative 
estimate for the validity of a hydrodynamic description is given by the condition that 
relaxation time is small compared to the time scale 1/u defined by average trap 
frequency w [H |9]. For contact interaction this condition can be expressed in terms 
of the scattering length to measure the strength of atomic interaction [HT]. The 
strength of dipolar interaction is usually measured in terms of the convenient length 
Ld = 3add = fJ'Op''^^/ (4vr^^), where fx is the magnetic moment of an atom of mass m and 
Ho the permeability of free space. Using this length scale to measure dipolar interaction, 
the condition for the validity of hydrodynamic description can be expressed as [9l [31] 



where N is the number of atoms, Oho = yh/{mui) the oscillator length, T the 
temperature, Tp the Fermi temperature and F{x) an universal function and can be 
taken of the order of unity [21] under experimental conditions. Hence for sufficently 
large N and/or sufficiently small a^o the system should enter the hydrodynamic regime. 
For ^^^Dy atoms of magnetic moment 10 fiB with fiB the Bohr magneton, the length 
Ld = 3add ~ 3 X 130ao where qq is the Bohr radius and for an oscillator length Oho = 0.25 
/im and N = 5000 the system is already in the hydrodynamic regime. For polar Fermi 
molecule ^°K^^Rb of electric dipole moment d = 0.566 Debye in the singlet rovibrational 
ground state [23], the length Ld = md'^ /{Airh'^eo) ^ 3 x 2000ao [1], where eo is the 
permitivity of free space, and the system should enter the hydrodynamic regime with 
larger Oho = 1 /^m (weaker trap) and much smaller N of few hundreds. 

At sufficiently low temperature the normal Fermi gas enters the degenerate phase. 
Most macroscopic properties, like shape, density, chemical potential, sound propagation 
etc., of this system can be described by the Landau hydrodynamical equations [27] . 
A velocity field is introduced as the gradient of a velocity potential $ of fiow, usually 
related to phase of the order parameter, by v(r, t) = V$(r, t), subject to the irrotational 
condition V x v(r, t) = 0. The continuity and the fiow equations are then given by [2Z] 

^''^''^^+V-Kr,t)v(r,t)] = 0, (2) 



dt 

: 0, (3) 



dv{r,t) ^ 



+ H{n,r) + V{r) 
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where n{r, t) is density at space point r and time t, ^(r) is an external trap, usually 
taken as 

V^(r) = ^mu;2^z/V + AV), (4) 

with uj the trap frequency and v and A are anisotropy parameters. The bulk chemical 
potential yu(n, r) is determined by the equation of state of the uniform Fermi gas of 
dipolar atoms and is given by 

fi{n, r) = + J dr'n{r')Udd{r - r'), (5) 

where the first term on the right-hand side (rhs) is the Fermi energy Ep = 
/i^(67r^n)^/^/(2m) [21] and the last term is the contribution from dipolar interaction 
energy [28], with f/(id(i" — r') is the dipolar potential. 

The hydrodynamic description is valid for a macroscopic observable with its 
characteristic excitation wave-length A much larger than the healing length [23]. A 
safe condition to satisfy this criterion is to take the wave length to be much larger than 
de Broglie wavelength at the Fermi surface, i.e., [231 132] 

A > 2n/kp (6) 

with Fermi momentum kp defined by Ep = h'^kp/2m. Truly speaking, a degenerate 
Fermi gas may not be fully irrotational and allows rotational components in the velocity 
field, not allowed in the present hydrodynamic formulation. This fact should influence 
only the rotational properties [23] of the degenerate Fermi gas not considered in this 
paper. We also assume the absence of any velocity dependent frictional force. 

An approximate TF profile for density can now be obtained by setting velocity 
V = in g, when [23] 

/i(n,r) + y(r)=/iJF, (7) 

where fi^^ is the chemical potential of the trapped gas. When ^ is solved subject to 
the appropriate normalization condition, we obtain both the chemical potential /ij ^ and 
the density n{r). 

There is an equivalent description of the trapped degenerate Fermi gas in the LDA, 
based on the assumption that, locally the dipolar Fermi gas would behave like a uniform 
gas, so that the energy density can be written as the energy of the uniform system times 
the local density [23]- For the degenerate dipolar Fermi system, the classical energy 
density per particle is given by 



5 2m 



+ N^I dv'n{v)n{v')UM{v - r'), (8) 

where n(r) is the density of fermions normalized as / (irn(r) = 1. The first term on the 
rhs of ([s]) is the total kinetic energy of the spin polarized fermions filling all levels up to 
the Fermi sea, the second term is energy in the trap, and the third term describes the 
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dipolar interaction. A minimization of the classical energy ([8j) leads to the TF condition 
0- 

As the degenerate Fermi gas is a quantum system, a quantum pressure term when 
included in (|8| yields the following expression for net energy density 

fW = — |V.v/^r)l' + ^ci(r). (9) 



This gradient correction term |33j to the TF energy density (|8j) takes into account the 
additional kinetic energy due to spatial variation of density (near the surface). Such a 
surface term was first considered by von Weiszacker [311 ESI ES] in the description a large 
nuclei. Previous descriptions of a degenerate Fermi gas considered different coefficients 
in this term |25l [33] . The energy density ^ has successfully used in many problems of 
Fermi gas [53 [351 [23 [M] • 

With the Lagrangian density C{r) = S{r) — Honijc) the Euler-Lagrange equation is 
given by [37] 



8m ^ ' 2m^ ^ 



+ nJ Uddiv - v')n{v')dv' 



n(r), 



(10) 



with /io the chemical potential. The derivarive term in (10) term contributes much 



less than the dominant "Fermi energy" term ^^[67r^A^n(r)]^/'^/(2m) and its neglect 
leads to the TF relation ([7]). The dipolar interaction in ( [Io[ ) is taken as Udd{^) = 
3arfrf^^(l — 3 cos^ 9) / (mR^), R = r — r' where 9 is the angle between the R and the 
polarization direction z. 

The condition ([T| refers to the validity of a hydrodynamic description of the system 



(v 7^ 0) P[3l]. A reliable stationary description obtained using (10) for v = 0, of density 
and other macroscopic properties, such as sound velocity, as considered in this paper, 
can be obtained for a smaller number of fermions consistent with condition ^ |2ll [32] , 
provided the contribution of the kinetic energy term ^^[67r^A^n(r)]^/^/(2m) in (10) is 
much larger than that of other terms. 

It is convenient to write a dimensionless form of ( 10 ) with the potential ^ as 



+ 3addN 



1 - 3 cos2 9 



i?3 



n{r')dr' 



n(r] 



:iii 



where energy, length and density are expressed in units of hoj, Iq = Jfi/muj and Iq . 



2.2. Variational Approximation 



The energy density corresponding to the dimensionless equation (11) can be written as 

SaddN f / /N 1 - 3 cos^ 9 

ar n{r)n{r 



Sir) 



1. 



IV. 



n(r] 



+ i(^y + AV)n(r) 



10 



(12) 
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A variational approximation for the problem can be obtained with the following 
Gaussian ansatz for density [30] 



n(r) 



exp 



wt 



(13) 



where Wp and Wz are the variational widths along radial p and axial z directions. With 
this density, the effective energy per particle of the system E = J drS{r) is 



E 



1 1 



w„ 



2wi 



+ 



+ 



Naddfjn) 

p ^ 



+ 



2 4 
/3 9(67r^iV)^/3 
5 bOmut^^wl^^ 



(14) 



where k = Wp/wz and 



1 + 3/s;2tanh"Vl - 



^2 



(1 - ^2)3/2 



(15) 



The variational equations are obtained by minimizing the energy (14) by dE/dw 
dE/dwp = HEO]: 

2 1 add Ng{K) , [3 6{6n^Nf/^ 



where 



9W 



h{K) 



4^3 



4w? 



27r w'^pW, 



2vr w'j^wl 




add 2Nh{K) 3 6(67r2iV) 



5 2hT^w'J'^Wp'^ 



2/3 



5 257rwf ^^if 



(16) 
(17) 



2 - 7/t2 - 4k'^ 9K^tanh-Vl - fi;2 



:i - K 



2^2 



+ 



(1 - ^2)5/2 



1 + 10k2 - 9K2tanh-Vl - 



(1 - /«2)2 (1 _ fi;2)5/2 

The chemical potential yUo of the system per particle is 



(19) 




2 Naddf{K 
2nWpWz 



+ 



2 4 

5 9(67r2Ar)2/3 

3 SOTrw^^Wz^'^ 



(20) 



2. 3. Approximate density for cigar and disk shapes 

In many situations of experimental interest, the Fermi gas could be subject to a strong 
trap either in the polarization z direction or in the transverse radial p plane. The 
system then has a one-dimensional (ID) cigar or 2D disk shape, respectively. In such 
cases simplified equations in lower dimensions could be very useful [39] . 

First we consider the reduced ID equation for a cigar shape. We assume that 
there is a strong trap in the x — y plane and that the density in this plane is given by 
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the Gaussian ground state [39j n{p) = exp{—p'^/dp)/{y/7idp),dp = ^J\/{2u) of the trap 
so that the 3D density n{Y) satisfies 



n[r) 



2dl 



(21) 



Substituting this density in (11), and multiplying this equation by the corresponding 



Gaussian yn{p) and integrating out the p dependence we obtain the reduced ID 
equation |Ml IH] 



PiD\/n{z) 



+ 



d^ 



dk 



sioiC) 

where n{kz) = e 



CO 2^ 

du- 



n{kz)s 



ID 



lOd'/' 
kzdp \ 



V2 



n{z), 



u 



(22) 
(23) 



M + 

'^^^^n{z)dz^ and piD is the chemical potential. An approximate 
variational solution of (22) is possible with the following ansatz for density [H] 
n{z) = l/{^/^lWz)exp[—z'^/wl], while the width is determined by solving (17) with 



dp and k = dp/wz 



Next we consider the reduced 2D equation suitable for a disk shape with a 
strong axial trap. The dipolar Fermi gas is assumed to be in the ground state [31] 
n{z) = exp{—z'^/dl)/{y/7idz), dz = ^JlJ{2X), of the axial trap X^z"^ /2 and the 3D density 
can be approximated as 



n{p) 



'/\/d. 



■ exp 



z 
'2rf! 



(24) 



Using this ansatz in (11), and multiplying by the corresponding Gaussian \Jn{z) and 
integrating out the z dependence we get the effective 2D equation [HI WI\ 



P2D^/n{p) - 



^ p 



2 2 I 

-z/ p + 



[3[6Nn{p)]^/^TT 
5 W 



d^k 



(2vr) 



n{kp)h 



2D 



kpdz 



n{p), 



(25) 



where n(kp) = / e^^fPn{p)d'^p,h2D{0 = 2 — 3^/7T^e^'^[l — erf(,^)]. An approximate 
variational solution of (25) is possible with the following ansatz for density [H] 
n{p) = 1/ {nw^) exp[—p'^w^], while the width Wp is determined by solving (16) with 
Wz = dz and k = Wp/dz- 



2.4- Sound propagation in a uniform Fermi gas 

To find the sound velocity in a uniform Fermi gas, we evaluate the Bogoliubov spectrum 
using the linearized hydrodynamic equations. We consider a stationary Fermi gas in a 
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box with periodic boundary condition with the trap V fixing just the allowed plane wave 
solution. Then (|3]), after the inclusion of the gradient term of reduces to 

f=i^ + Kn, r) + m^^^ = 0. (26) 

Now we allow small perturbation in n and $ around their equilibrium values uq and $o 
by n(r,t) ^ tiq + n(r,t) and $(r,t) ^ $o + $(r,t), then ^ leads to 

c?n(r, t) 



+noV'<l>(r,t) =0. 

In (26|, we need to use ^Jn{r, t) ^ ^/uq + n{v,t)/{2^/n^) and [n(r, t)]^/^ 
2n(iv?)/(3no"^'''^), while /i(n, r) ^ jlo + fi{n,r) with 

;i2(67r2^2/3, 



n, 



(27) 

2/3 _^ 



/i(n,r) 



+ J dr'n{r',t)Uddir - r') 



where fxo is the stationary value of fx. Then ( [26| leads to 
+^ n,r +m ^ ' = 0. 

Assuming the perturbations n and $ in plane-wave forms n = f;,o6xp[z(k • r) 
$ = $0 exp[i(k ■ r) — wt], /i has the form /2 = no/xo exp[i(k ■ r) — ut] with 



;i2(67r 



2\2/3 



3mn, 



1/3 



+ 



m 



(3 cos^ 



(28) 

(29) 

ut] and 

(30) 



where ^ is the angle between the propagation direction and polarization direction z, and 
the last term in (30) is just the Fourier transform of the dipolar potential. Then (27) 
and (29) become: 



iujriQ — noA;^$o = 0, 



no — iojm^Q = 0. 



(31) 
(32) 



16mno 

The condition of existence of the solution to this set of equations leads to the Bogoliubov 
spectrum 



efc = hu 



^ Am 



Am 



+ Ariofio 



(33) 



The sound velocity, defined as Cs{9) = limk^Q^ek/hk) for a uniform dipolar Fermi gas 
can be written as 



+ 



AixriQaddh'^ 



(3cos2^- 1) 



(34) 



where vp = hkF/m = h^GTr'^noY^^ /m is the Fermi velocity. For a nondipolar Fermi gas 
{o-dd = 0), (34) leads to the well-known velocity of vp/y/S 



The angle-dependent second term under the square root in (34) is responsible for 
anisotropic sound velocity. Specifically, for 9 > 54.73° degrees, this term is negative 
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corresponding to a decrease in sound velocity. For large dipolar interaction, for 6 greater 
than a critical value and for large density no the sound velocity could be imaginary 
corresponding to no propagation. However, in this study we shall only consider moderate 
values of density and dipolar interaction, that allow anisotropic sound propagation in 
all directions. 



3. Numerical Results 



3.1. Stationary properties of trapped dipolar Fermi gas 



For a trapped 3D Fermi system we solve (11) numerically after discretization The 



divergence of the dipolar term at short distances has been handled by treating this term 



in momentum (k) space |6]. For numerical calculation in section 3.1, we consider a 
degenerate Fermi gas of ^^^Dy atoms with add ~ 130ao and employ the oscillator length 
/o = 1 /xm. 

The anisotropic dipolar interaction is partially attractive in certain angles and 
repulsive in others and contributes very little in a spherically symmetric trap. The 
dipolar interaction contributes attractively in the cigar-shaped configuration along the 
polarization z direction and repulsively in the disk-shaped configuration perpendicular 
to the polarization z direction. Hence we will mostly consider the degenerate dipolar 



Fermi gas in cigar and disk shapes. The 3D model (11) is very strongly nonlinear with 
nonlinearity J\f = (67r^A^)^/^/2, leading to large nonlinearities oi J\f = 760 and 3526 for 
= 1000 and 10000, respectively. 
First, we compare the results of the reduced 2D density n{x,y) of a disk-shaped 
degenerate Fermi gas of ^^^Dy atoms with A = 10, = 1 as obtained from the 3D and 



2D models (11) and (25), respectively, and from the variational approximation to the 3D 
model. The reduced 2D density n{x,y) is defined by n{p) = n{x, y) = dzn{x, y, z). 
In figure [l] (a) we plot the reduced 2D density along x axis n(x,y = 0) for different N. 
Next, we compare the results of reduced ID density n{z) of a cigar-shaped degenerate 
Fermi gas of ^^^Dy atoms with A = 0.1, u = 1 as obtained from the 3D and ID models 



(11) and (22), respectively, and from the variational approximation to the 3D model. 
The reduced ID density n{z) is defined by n{z) = J^^dx J^^dyn{x,y, z). In figure [T] 
(b) we plot the reduced ID density along z axis n{z) for different A^. From figure [l] we 
find that in both cigar and disk shapes the ID and 2D models perform fairly well, even 
for large nonlinearities, when compared with the results of the full 3D model. To see if 
the condition ^ is satisfied by the shapes in figure [l| we can use the TF estimate for 
kp of a trapped degenerate Fermi gas [24|: kp ^ {A8NY^^/Iq. For 100 < < 10000, as 
in figure [l| the de Broglie wave length at the Fermi surface is 271 /kp Iq ■ The shapes 
and sizes in figure [l] are much larger than this value, consistent with the condition ([6]). 

Next we consider a cigar-shaped degenerate Fermi system of ^^^Dy atoms with trap 
parameters = 1, A = 0.3. In figure [2] (a) we plot the numerical and variational results 



for the chemical potential yU and rms sizes {x) and {z) of this system for 10 < < 10000. 
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Figure 1. (a) Radial density of a disk-shaped Fermi gas of ^^^Dy atoms in x ~ y 
plane n{x,y = 0) = f^^dzn{x,y = 0,z) from the 3D (111 for = 1,A = 10 
compared with its variational (var) approximation and the solution of the reduced 
2D ([25| for iV = 100, 500, 10000. (b) Axial density of a cigar-shaped Fermi gas of 
^^^Dy atoms along polarization z axis n{z) = I^oo dxdyn{x, y, z) from the 3D 
(111 for u = \,\ = 0.1 compared with its variational (var) approximation and the 
solution of the reduced ID ^ ioi N = 100, 1000. 



In figure [2] (b) we plot the same for a disk-sfiaped degenerate Fermi system witli trap 
parameters = 1,A = 10. Finally, in figure [2] (c) we plot the same quantities for 
= 1000 versus the trap anisotropy A for v = 1. In all cases the variational results are 
in good agreement with the 3D model. For medium to small number {N < 10000) of 
trapped ^^^Dy atoms as considered here and also of experimental interest, the effect of 



the dipolar term in Eq. (11) (the last term in this equation) is small compared to the 
Fermi energy term (the last but one term there). Hence the effect of the dipolar term 
in figures [2] is small. The plots in these figures only change by less than about two to 
four percent if we set a^d = 0. 

3.2. Anisotropic sound propagation in uniform dipolar Fermi gas 



The hydrodynamic analytical result for sound velocity (34) in a uniform dipolar Fermi 



gas shows a clear anisotropy through the angle 6 between the propagation direction 
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Figure 2. (a) Reduced chemical potential fi/{'huj), and rms sizes (x), and (z) of a 
cigar-shaped Fermi gas of ^^^Dy atoms with = 1, A = 0.3 versus N from a numerical 
(n) solution of the 3D (11) and variational (v) approximation [ (161, (17) and (20)]. 
(b) The same for a disk-shaped gas with ly — 1, X — 10. (c) The same for the trapped 
Fermi gas of 1000 ^^^Dy atoms versus the parameter A with v — \ from a numerical 
solution of (11) and variational approximation. 



r and the polarization direction z. In this expression for sound velocity there are two 
competing terms under the square root. The first term ff'/S involving the Fermi velocity 

2/3 

is isotropic and proportional to Uq whereas the second term proportional to dipolar 
interaction is anisotropic and proportional to density no- The anisotropy in sound 
velocity Cs{0) will manifest strongly for large strength add and for large as the 
anisotropic dipolar term in (34) will increase more rapidly with no than the isotropic 



term. In figure [3] (a) we plot the sound velocity Cs{9) for angles 6' = and vr as a 
function of density tlq for a degenerate Fermi gas with add = 130ao (^^^Dy atom) and 
2000ao (polar "^^K-^^Rb molecules in the singlet rovibrational ground state P) using the 
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analytical result (34). The sound velocity for nondipolar atoms {add = 0) is also plotted. 
From this plot we find that the anisotropy in sound propagation as measured by the 
difference [cs(0) — Cs(7r/2)] is sizable only for a relatively large density no = 10^^ cm~^ 
for ^^^Dy atoms. However, for polar ^"K-^^Rb molecules the anisotropy is appreciable 
for a relatively low density of Uq = 10^^ cm~^. From figure (a) we see that for ^''K-^^Rb 



molecules, the sound velocity is imaginary for 6' = for a density of about no = 5 x 10^^ 



cm ^. 



To study sound propagation, the present stationary (static) 3D model ( pAj ) is 
generalized to include time variation by replacing /i by the usual time derivative ihd/dt. 





10 20 30 40 50 60 70 80 

6 (degree) 



Figure 3. (a) The analytical sound velocity Cs{9) given by (34) in axial z (9 = 0) 
and radial p {6 — 7r/2) directions versus density for add = 130ao and 2000ao as 
well as the velocity for add ~ 0. (b) The analytical (ana) sound velocity Cs{0) in axial 
z {0 = Q) and radial p (9 = 7r/2) directions versus add compared with the results of 
numerical simulation (num). (c) The analytical sound velocity Cs{9) versus polar angle 
9 compared with the results of numerical simulation. 
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Consequently, the infinite dipolar Fermi gas satisfies the following Galilei-invariant 
equation [37] 



ih 



d 



dt 
+ I2a. 



2(2m) ' 
1 - 3 cos^ 9 



2V{r) 



77-0 (r') fir' 



2|-[6.X(r)]^/3 



^o(r), 



(35) 



2m J i?3 

where density no(r) = Nn{r) is not normalizable for infinite hydrodynamics. Equation 
(35) is consistent |37] with the time-dependent hydrodynamic equation ^ after the 
inclusion of the gradient correction term [33] . 

The anisotropy of the dipolar interaction would be prominent at low to medium 
density for large dipolar interaction. To illustrate the anisotropic sound propagation 
we will consider two examples: ^^^Dy atoms at a medium density of 10^^ cm~^ and 
the polar molecules ^°K-^^Rb at the low density of 10^^ cm~^. First we consider the 
numerical simulation of sound propagation in the infinite ^^^Dy gas {add = 130ao) at 
a background density of uq = 10^ /im"^ = 10^^ cm~^. The numerical simulation is 
initiated with a 3D Gaussian pulse at the center of the uniform 3D background density 
given by no{r) = (10'^ + lO^e"''^/^"'^) /im~^, w = 2 /im, subject to a weak expulsive 
Gaussian potential ^(r) = O.OOOOle"''^/^"'^ /xm"^. With this initial condition (35) is 



16 i 
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Figure 4. The contour plot of 3D density in x ~ z plane nQ{x,y = 0,z) during 
sound propagation on a background density no = 10"^ ^m^^ of a degenerate dipolar 
i^^Dy gas with add = 130ao at times t = (a) O.lSio, (b) 0.2io, (c) 0.25io, and (d) 0.3^0 
(to = 0.005 s). The contour plot of 3D density in x — y plane nQ{x,y,z = 0) during 
the same sound propagation at times (e) 0.15to, (f) 0.2to, (g) 0.25to, and (h) 0.3io- 



Study of a degenerate dipolar Fermi gas of ^^^Dy atoms 



15 



to 

N 

-8 

-15 



16 
_ 8 

to 





® 


® 




(a) ' ' ' (b] ■ ■ (c) ■ ' ' (d) ■ 

-8 816 -8 816 -8 816 -8 816 
X (nm) X (^im) X ((im) x (^im) 


• 


® 







1 (nm"'^) 
9.95 (^m"^) 



1 (tim"'^) 
9.95 (>im"^) 



-8 816 -8 815 -8 816 -8 816 
X (nm) X (^im) X ((im) x (^.m) 



Figure 5. The contour plot of 3D density in x — z plane no{x,y = Q,z) during 
sound propagation on a background density uq ~ 10 ^iycT^ of a degenerate dipolar 
molecular "OK-^^Rb gas with = lOOOoo at times t = (a) 0.2io, (b) 0.5io, (c) O.Sio, 
and (d) l.lto (io = 0.004 s, Zq = 1 ^m). The contour plot of 3D density m.x~y plane 
no{x,y,z = 0) during the same sound propagation at times (e) 0.2to, (f) 0.5to, (g) 
0.8<o, and (h) I.I^q. 



solved by real-time propagation [43j to study sound waves. An ellipsoid-like sound wave 
front is found to emerge outwards upon time propagation. From a numerical study of 
this wave front the sound velocity in different directions is calculated. 

Typical anisotropic sound propagation in a ^^^Dy degenerate gas of density 10^^ 
cm~^ is shown in figure |4j The sound propagation is illustrated via contour plots of 3D 
density nQ{x,y,0) in the x — y plane and of density no{x,0,z) in the x — z plane. In 
figures |4] (a), (b), (c), and (d), to illustrate the anisotropic sound propagation in the 
X — z plane, we show the contour plot of nQ{x,0,z) at times t = 0.15to, 0.2to, 0.25to 
and 0.3to, respectively, with to = 2mll/h'^ ~ 0.005 s the time scale and Iq = 1 /im the 
length scale used in the numerical solution of (35). However, the propagation in the 
radial x — y plane is isotropic. This isotropic sound propagation is shown in figures |4] 
(e), (f), (g), and (h) via the contour plot of no{x,y,0) at times t = 0.15to, 0.2to, 0.25to 
and 0.3to, respectively. A clean wave front of high density is identified in these contour 
plots encompassing a region of low density — a typical panorama in sound propagation. 

The sound propagation along z direction (polar angle ^ = 0) has a velocity 
(cs = Vp/y/S) larger than that for a nondipolar system {add = 0) of same density, 
whereas that in the x — y plane (polar angle 6 = 7r/2) has a velocity smaller than that 
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for a nondipolar system as can be seen from (34). The analytical sound velocity Cg for 
nondipolar atoms of density 10^^ cm~^ and atomic mass 161 is 8.88 mm/s compared to 
the numerically obtained velocity of 8.5 mm/s. For ^^^Dy atoms of density no = 10^^ 
cm~^ with add = 130ao, the analytical radial sound velocity is Cs(7r/2) = 8.1 mm/s 
(numerical 7.0 mm/s), and the analytical axial sound velocity is Cs{0) = 10.3 mm/s 
(numerical 10.0 mm/s). In figure [3] (b), we present the variation of axial {9 = 0) and 
radial {6 = 7r/2) sound velocities versus add as obtained from numerical simulation 
and analytical consideration (34). In figure |3] (c), the analytical result for velocity 
versus the polar angle 6 is compared with the numerical result. The agreement between 
the analytical result (34) for sound velocity and the result of numerical simulation is 
satisfactory considering the very large nonlinearities present in the system due to large 
density (no) and large dipolar interaction. 

To have a larger anisotropy in sound propagation than in a gas of ^^^Dy atoms 
at a lower density a stronger dipolar interaction, as in polar molecules, is needed. A 
convenient polar fermionic molecule molecule ^°K-^^Rb has a permanent electric dipole 
moment 0.052 Debye {add ~ 20ao) for the triplet rovibrational ground state and 0.566 
Debye {add ~ 2000ao) for the singlet rovibrational ground state [1]. For illustration, 
we consider a uniform gas of ^''K-^^Rb molecules of density hq = 10^^ cm~^ with 
O'dd = IOOOoq. The numerical simulation is initiated with a 3D Gaussian pulse at the 
center of the uniform 3D background density given by no(r) = (10 + e-''^^^'"^) /im-^ 
w = 2 yum, subject to a weak expulsive Gaussian potential V{r) = O.OOOOle"''^/^"'^ 
/zm"^. With this initial condition (35) is solved by real-time propagation [33] to study 
sound waves. Typical anisotropic sound propagation in a ^°K-®^Rb degenerate gas of 
density 10^^ cm~^ is shown in figure 5. Because of the stronger dipolar interaction, a 
larger anisotropy has appeared at a lower density in figure 5 compared to figure |4| For 
^°K-^^Rb molecules of density uq = 10^^ cm~'^ with add = IOOOoq, the analytical radial 
sound velocity is Cs{tt /2) = 2.05 mm/s (numerical 2.4 mm/s), and the analytical axial 
sound velocity is Cs{0) = 3.03 mm/s (numerical 3.4 mm/s). 



4. Summary 

We developed a 3D theoretical formulation for describing certain macroscopic 
observables of a degenerate dipolar Fermi gas appropriate for studying stationary 
properties, such as, shape, size, energy, chemical potential etc. of a trapped system. 
The effect of dipolar interaction is negligible in the spherically-symmetric configuration 
and the dipolar interaction manifests strongly in the asymmetric cigar and disk shapes. 
Simple reduced equations in ID and 2D suitable for studying the trapped degenerate 
dipolar Fermi gas in cigar and disk shapes, respectively, are derived. Also, a Gaussian 
variational approximation for studying these macroscopic properties is derived. We 
apply the present formalism to study the stationary properties of a trapped degenerate 
Fermi gas of ^^^Dy atoms. The stationary properties of the 3D model under appropriate 
conditions are found to be in satisfactory agreement with those from the reduced ID 
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and 2D models as well as with variational approximation. 

Using the present 3D model we also obtained analytical results for anisotropic sound 
velocity as a consequence of the anisotropic dipolar interaction in a dipolar Fermi gas in 
agreement with numerical simulation of the 3D (35). The sound velocity is larger along 
the polarization direction than in the transverse plane. 
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